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We consider models of quasi- 1-d, planar atomic wires consisting of several, laterally coupled rows 
of atoms, with mutually non-interacting electrons. This electronic wire system is coupled to phonons, 
corresponding, e.g., to some substrate. We aim at computing diffusion coefficients in dependence on 
the wire widths and the lateral coupling. To this end we firstly construct a numerically manageable 
linear collision term for the dynamics of the electronic occupation numbers by following a certain 
projection operator approach. By means of this collision term we set up a linear Boltzmann equation. 
A formula for extracting diffusion coefficients from such Boltzmann equations is given. We find in 
the regime of a few atomic rows and intermediate lateral coupling a significant and non-trivial 
dependence of the diffusion coefficient on both, the width and the lateral coupling. These results, 
in principle, suggest the possible applicability of such atomic wires as electronic devices, such as, 
e.g., switches. 
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I. INTRODUCTION 

The strong desire for an ongoing miniaturization of 
electronic circuits has led to the idea of molecular elec- 
tronics. Very roughly molecular electronics include all 
implementations of basic electronic functional devices on 
the scale of atoms or molecules. The most basic device 
surely is a wire, meaning a system in which an electronic 
current may flow as a response to an electric field in a 
prescribed direction. Various implementations of such 
wires have been suggested and experimentally explored. 

Some suggestions are based on chemically synthesized, 
mechanically more or less stable chainlike molecules or 
atomic structures [H-ll, @| j which are contacted to some 
kind of leads at both ends, but otherwise " free- hanging" 
without, e.g., connection to some substrate. There is 
experimental data available on the conductivities of wires 
consisting of a single molecule 0, @. So far molecular 
wires of lengths up to several nm could be experimentally 
produced. 

Furthermore, there are various suggestions and realiza- 
tions of wire-like, 1-d, nanostructures based on carbon, 
e.g., carbon nanotubes or graphene-nano ribbons, placed 
on some substrate. Those carbon-systems show remark- 
able electronic and transport properties (see, e.g., [llj 
and references therein) . Both of them may feature either 
metallic, semimetallic or semiconducting band structures 
depending on chirality and width, respectively. 

Some other experimental of quasi- 1-d wires of nm- 
scale-diameter feature crystalline cores. Chemically sta- 
ble silicon or metallic nanowires of that type have been 
successfully fabricated d, [l3 ■ It has also been shown by 
experiments Q (scanning tunneling spectroscopy) and by 
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theory (density functional theory) Q that the electronic 
structure of those nanowires is similar to the electronic 
structure of the corresponding bulk material for diame- 
ters larger than cx 7 nm. However electronic structure 
and resistivity may differ from bulk-values for smaller 
diameters. 

Yet another scheme envisions the formation of chains 
or rows of conducting atoms or molecules on the surface 
of some insulating substrate. The formation of such lin- 
ear molecular structures is thought to be achieved by self- 
assembly, possibly influenced or controlled by structures 
on the substrate surface. First steps in that direction 
have been demonstrated [12 1. 

Two points should be stressed here: i. So far none 
of the above techniques provides atomic wires which are 
scalable, stable and sufficiently controlled to set up com- 
plex electronic circuits, ii. Assume, nevertheless, arbi- 
trarily long, stable wires were obtainable, then, accord- 
ing to all the above suggestions, except for the first "free 
hanging" one, the electronic system would be (weakly) 
coupled to bulk-type phonons. The latter being either 
due to the lattice of the crystalline wire itself or simply 
to the substrate. Even if impurities and electron-electron 
coupling were negligible the phonon-coupling would in- 
evitably control the transport properties. 

In the work at hand we essentially aim at analyzing 
the peculiarities of the transport properties of some at 
present hypothetical, long atomic wires, the electrons of 
which are subject to scattering along the wire as men- 
tioned under the above second point. 

Theoretically electronic transport through atomic 
wires is often analyzed by using a Landauer-Biittiker- 
approach 7]. Basically, the latter addresses a scenario of 
electrons moving ballistically, i.e., unscattered, through 
an ideal 1-d conductor from one electrode to another (due 
to an applied voltage) . The resulting (finite) conductivity 
is mainly determined by the properties of the contacts. 
While this method is well suited for the above, relatively 



2 



short, "free hanging" wires the incorporation of rather 
complex scattering along the wire appears challenging. 

A standard approach to transport in bulk electron- 
phonon coupled systems is the Holstein approach. Here 
one relies on an evaluation of a Kubo-formula by means 
of an approximate diagonalization of the system based on 
a transformation into a polaron picture. However, this 
produces, in the limit of weak coupling relevant, diverg- 
ing terms which would render the coupled model ballistic. 
This can only be cured by the ad-hoc introduction of a 
phenomenological (isotropic) polaron lifetime, the latter 
being due to either impurity or electron-electron scat- 
tering, i.e., scattering mechanisms which are not present 
in the pure electron-phonon model [13]. Altogether this 
approach appears inadequate for the analysis of weakly 
coupled but highly anisotropic electronic systems like the 
atomic wires addressed below. 

Our atomic wires may be viewed as hypothetical, ideal 
realizations of the "self-assembled-atoms- type" . Our in- 
vestigations are based on substantially simplified models. 
We consider structures consisting of several, parallel, in- 
finite mono-atomic chains of atoms. Loosely speaking 
the electrons may hop between the atomic sites, possibly 
with different hopping strengths for the lateral and the 
longitudinal direction. The band structure of such a sys- 
tem may be of the metallic, insulating or semi-conducting 
type, (just like the graphene nano-ribbons), depending 
on the number of rows and the hopping strengths. We 
take the electrons as mutually non-interacting but weakly 
coupled to some bulk phononic system. The sorround- 
ing system is assumed to be insulating, i.e., there is no 
electronic transport from the wire into the sorrounding 
system. (While non-interacting electrons are surely an 
idealization, this idealization may be adequate for 1-d 
structures which tend to remain ballistic, even in the 
presence of interactions [l3 - [T7| ). Since, as explained 
above, neither the Landauer-Biittiker nor the Holstein 
approach appear well-suited for such a system we resort 
to a yet different approach based on Boltzmann equation 
descriptions. 

This paper is organized as follows. In Sec. we first 
introduce the above quantum model of the wire in de- 
tail. In Sec. mil we formulate, based on the model, a 
linear (ized) master equation for the dynamics of some 
variables corresponding to occupation numbers of elec- 
tronic momentum modes. This is done on the basis of an 
approach derived in detail in [ll]. In Sec.|TV]we suggest 
a Chapman-Enskog-type calculation of a diffusion coeffi- 
cient from a linear Boltzmann equation. By means of the 
previously formulated master equation we set up such a 
linear Boltzmann equation describing the electronic dy- 
namics of our model. Finally, in Sec. |V] we numerically 
evaluate the obtained formula for the diffusion coefficient 
of our model. We thus analyze the dependence of the dif- 
fusion coefficient on the number of neighboring rows as 
well as on the lateral coupling strength between neigh- 
boring rows, finding strongly discontinuous behavior. We 
close with summary and outlook. 



II. INTRODUCTION OF THE MODEL 

Routinely we open the considerations by specifying 
the underlying quantum model on the basis of which we 
below analyze the dynamical properties of the atomic 
wires. The Hamiltonian consists of three parts: an elec- 
tronic part -ffo.oi describing the hypothetically unper- 
turbed electrons on the wire, a phononic part ifo.ph rep- 
resenting a standard set of decoupled phononic modes 
and an interaction V which couples the electrons to the 
phonons: 

B 

} r=l i 

-ffo.ol -ffo.ph 

+ (E E ^-Ur-K.K + h.C.), (1) 

k,q r.s=l ^ 
^ ^ ^ 

V 

with n = BL. 

A real space sketch of the electronic part is depicted 
in Fig.[TJ We assume a quasi-l-d, planar, cubic lattice of 
atoms that forms the wire, L corresponds to the length. 
Concretely we assume that the electrons may occupy lo- 
calized (Wannier) states at those lattice sites, i.e., a stan- 
dard tight-binding description. The wire now consists of 
a finite, possibly small integer number B of parallel rows 
of atoms. The motion of the electrons along the rows is 
controlled by nearest neighbor hopping matrix elements 
of strength Tn . The hopping from one row to the adja- 
cent row(s) is controlled by matrix elements of strengths 
T±. The electrons on the wire are assumed to experi- 
ence no mutual interactions, i.e., without any further cou- 
pling their dynamics would be simply ballistic. Of course, 
the tight-binding-description for the non-interacting elec- 
tronic system is a strong idealization, which, however, 
may be a good approximation for 1-d systems. Since 
the wire is furthermore strictly periodic along the rows, 
its Hamiltonian may be denoted based on the dispersion 
relation e^^r as done in ([ij. Here ej^^ denotes the one- 
particle-energy corresponding to the lattice-momentum j 
and the "branch" r. Due to the B parallel rows the dis- 
persion relation features B branches or "bands" . Fourier 
transforming the Hamiltonian described by Fig. [T] essen- 
tially yields 

e^^r = -2T|| cos j + T^Er , (2) 

where the concrete have to be determined from the 
subsequent diagonalization of a BxB-Matrix. (For an 
B — 2 example see Fig. [21) Hence, the bandwidth of 
each band (which we label Eb) is determined by the lon- 
gitudinal coupling Eb = 4T||, whereas the displacement 
T±er scales linearly with the lateral coupling T±. That 
is, the model may feature a metallic, semiconducting or 
insulating band structure, depending on the choice of the 
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system parameters, namely the hopping strengths T\\ , Tl. 
As it is demonstrated below, increasing Tl for fixed T\\ 
allows for a transition from metallic over semiconducting 
to insulating behavior. 

The phononic dispersion relation cji may take some 
standard form for acoustic or optical phonons. Below it 
will turn out that the details are irrelevant within the 
frame of our considerations as long as phononic energies 
are much smaller than electronic energies, i.e., Wi << ej^r- 
The latter will be assumed. 

The interaction V is also of standard lowest order form, 
it may be viewed as corresponding to an electron scatter- 
ing under creation or annihilation of a phonon, with the 
phonon essentially compensating for the electronic mo- 
mentum change. In our approach the electron-phonon- 
coupling V is treated as a weak perturbation to the non- 
interacting electronic system. 
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FIG. 1: Sketch of the quantum model we use to describe the 
properties of an atomic wire. Circles denote localized elec- 
tronic states on a regular (quasi-) 1-d lattice in configuration 
space. Lines indicate couplings with strength parallel and 
strength T± perpendicular to the wire. 



III. FORMULATION OF A MASTER 
EQUATION FOR COARSE GRAINED 
OCCUPATION NUMBERS 

As a first step in our analysis we now aim at mapping 
the quantum dynamics of some set of "coarse grained oc- 
cupation number deviations" in momentum space onto 
a master equation. To this end we follow directly the 
scheme described in [l8| . For an overview over mappings 
of quantum dynamics onto Boltzmann equations also see 
[Tsj and references therein. This scheme produces a linear 
master equation for a finite number of variables. These 
variables are sums of deviations of electronic occupation 
numbers from their equilibrium values. To make this well 
defined we have to specify which occupation number de- 
viations should contribute to a certain sum. Only with 
this specification the scheme from 18] may be imple- 
mented. We call all modes, the occupation number de- 
viations of which contribute to a certain sum, a "grain" . 
Thus the number of variables as well as the concrete rates 
from the above master equation depend on this graining. 
Since this graining is to some extent arbitrary, we must 
make sure that eventually our diffusion coefficient con- 
verges for a systematic refinement of the graining. This 



convergence then in turn sets the scale for a "correct" 
graining. In the following we first describe our choice 
for the graining and then give precise definitions for the 
above dynamical variables of the master equation. 

Our graining is realized by primarily partitioning the 
momentum space into M "energy shells" of equal en- 
ergy width /S.E. We then define the grains within one 
shell by means of the intersections of the electronic dis- 
persion relation with the respective energy shell as de- 
picted exemplarily in Fig. O We label the entity of all 
grains serially by Greek indices (e.g., k). (This may cor- 
respond to a consecutive numbering from left to right (see 
Fig. within one shell, which is then accordingly con- 
tinued through the shells with increasing energy.) Note 
that these grain indices already contain the information 
to which band the respective grain belongs, that is, the 
band indices r are absorbed into the grain indices. As 
mentioned above the necessary refinement of the graining 
is eventually determined by the convergence of the dif- 
fusion coefficient. However, we generally focus on grains 
that are small enough to allow for a linearization of the 
dispersion relation on each grain. 
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FIG. 2: Sketch of our specific coarse graining in momentum 
space for a model featuring B = 2 rows, i.e., the dispersion re- 
lation consists of 2 bands. Parallel dashed lines define energy 
shells of equal energy width AE. Rectangular boxes mark 
the individual grains. The Fermi edge is indicated by a thick 
line. Some shells contain grains from one band, some from 
both bands. Other model parameters: Eb = 3eV, T± — Ty , 
M = 10 energy shells and hah filling (e/ = 2.25eV). 

Now we give a mathematical definition for the above 
deviations of electronic occupation numbers from their 
equilibrium values which are summed over a grain to form 
a dynamic variable of the above master equation. An 
operator representing such a deviation is given by 



(3) 



where fj^r is the equilibrium value of the electronic occu- 
pation number. Since we are interested in weak electron- 
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phonon coupling, we may take fj,r from the Fermi distri- 
bution for the occupation numbers of the isolated elec- 
tronic system. Since we are furthermore interested in 
low temperature regimes, we may routinely replace the 
chemical potential /i in the corresponding Fermi distri- 
butions /j,r(A<,T) = (exp((ej,r - ^i)/kBT) + 1)^1 by the 
Fermi energy e p ■ With the help of such "individual mode 
deviation operators" Aj.,. the deviation operators for the 
whole grains A'' may simply be defined by the corre- 
sponding sums 



(4) 



where the sum of course runs over all modes belonging to 
some grain. The dynamical variables of the above mas- 
ter equation are then defined as the expectation values of 
the grain deviation operators d'^{t) for the actual quan- 
tum state p{t) of the full system evolving according to 
the coupled Hamiltonian H, i.e., d'^{t) := Tr{p(t)A''}. 
Naturally, the dynamics of p{t) are not directly accessi- 
ble. But, as explained above, the considerations in [l^ 
(which we, for brevity, do not intend to repeat or discuss 
here) provide a scheme by means of which the dynamics 
of the d'^ {t) may be mapped onto a linear master equa- 
tion. The scheme is believed to be valid as long as the 
electron-phonon coupling may be considered a small per- 
turbation to the electronic system Hq ^i- Thus, simply 
applying this scheme casts the dynamics into the form: 



(5) 



The concrete expressions of the rates may also be read 
off from the results of jl^ . As mentioned above we focus 
here on the case of negligible phonon energies which is 
explicitly treated in [18]. Exploiting this leaves us with 



dr 



■E 



|W^rs(k-i)|^ 



n 



(l+2gk-i) 



•cos[-(ei^s-ek,r)T], 

(6) 

where A''^ denotes the number of states within the grain 
K and 5q is the equilibrium occupation number of the 
phonon mode q, as taken from the Bose distribution 
for the isolated phononic system i/o,ph- (Note that the 
phononic occupation numbers are not treated as dynam- 
ical variables in our approach.) Basically, the rates ap- 
pear as integrals over some correlation functions, which 
usually decay within some correlation time t^, we gener- 
ally assume the latter to be short compared to the time 
scale of the relaxation dynamics (weak coupling limit). 
Further, we define 



r2 = 

riK. 



E 



|T4^,,(k-i)|2(l + 2.9k- 
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(7) 



Using this abbreviation and routinely exploiting the 
properties of the smc-function as well as the lineariz- 
ability of the dispersion relation on each grain, we can 
approximately perform the integration in ^ , finding for 
times larger than the correlation time for the rates 
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The 6E{ri},E{K.) is meant to indicate that transitions only 
occur between grains from the same energy shell, i.e., 
all other rates vanish. Note that the resulting rates are 
completely independent of the electronic equilibrium oc- 
cupation numbers, i.e., scattering above, below, or at 
the Fermi-level is in this approximation qualitatively the 
same. 



IV. CALCULATION OF A DIFFUSION 
COEFFICIENT FROM A LINEAR BOLTZMANN 
EQUATION 

Eventually we intend to calculate the diffusion coeffi- 
cient from a type of linear Boltzmann equation. To this 
end we are simply going to interpret the rates from (|5]) 
as the entries of the corresponding collision term. But 
even if the Boltzmann equation is fully defined, we have 
to specify how to calculate a diffusion coefficient from it. 
There are several approaches to the derivation of diffu- 
sion coefficients from linear(ized) Boltzmann equations 
p^l - [2^ . However, since we are here confronted with a 
subtlety that may not occur in other cases, namely a 
non-symmetric rate matrix, we suggest an alternative 
approach to the diffusion coefficient. This approach is 
outlined in the following. 

We start from a linear Boltzmann(-like) equation de- 
scribing the dynamics of some probability distribution 
Pni^, t) in a velocity discretized one-particle phase space: 

p,^{x,t) = -t;,,VxP,,(x,t) + Ei?,,„pK(x, i) . (9) 

The two addends represent the drift and the collision 
term, respectively. The velocity v appears as depending 
on some index rj. Eventually this index will be iden- 
tified with the quantum mechanical momentum mode. 
After Fourier transforming with respect to x we find 
for the (Fourier transformed) probability distribution 
Pv — Pr,(q,i) 



Pr, = y^^j-iqVijSjjK + R7jk)Pk , 



(10) 



where the drift term i?i now enters as a diagonal matrix 
and the collision term i? as a non-diagonal matrix with 
respect to the above representation. 

Eq. (fTU]) is formally identical to a non-specific, abstract 
time dependent Schrodinger equation where Ri + R cor- 
responds to some Hamilton operator. However, since 
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Ri + R is not Hermitian, eigenvalues may (and will) 
not be real. Furthermore, eigenstates represent spatial 
particle-density waves with wavevectors q which decay 
mono-exponentially. Such a behavior is in accord with 
diffusive transport. The corresponding decay constant is 
given by the imaginary part of the corresponding eigen- 
value. For arbitrarily small q (long wavelengths) one 
may treat the drift term Ri as a small perturbation to 
the collision term R. Therefore, in order to find the decay 
constant, we may calculate corrections to eigenvalues of 
R induced by i?i via a procedure that is essentially anal- 
ogous to standard time- independent perturbation theory 
as given in textbooks. In doing so, the eigenvectors of 
R represent the eigenstates of the unperturbed Hamilto- 
nian. The eigenvector Pq of R featuring the eigenvalue 
zero corresponds to the total, isotropic equilibrium dis- 
tribution in the above phase space. Thus, assuming that 
the evolving distribution that corresponds to a diffusion 
process exhibits locally (with respect to real space) more 
or less the equilibrium distribution with respect to veloc- 
ities (momenta), we have to analyze the decay dynamics 
of Pq induced by i?i . In other words we have to calculate 
the correction to the eigenvalue of Pq. Before going in 
any detail we note that the second order correction scales 
as (X q^. Thus, in case the first order vanishes, the decay 
constants for the corresponding density waves will scale 
as oc q^ as well. The latter dynamics represent a solution 
to a diffusion equation with a diffusion coefficient given 
by the prefactor of the q^-term of the decay constants. 
Thus, calculating the second order correction to Pq may 
and will turn out as being equivalent to calculating the 
diffusion coefficient. 

Nevertheless, as already mentioned above, we en- 
counter here the difficulty that the collision term R is 
in general non-symmetric and therefore the eigenvectors 
of R, denoted by {Pi}, are not necessarily orthonormal. 
This problem may be dealt with by introducing a dual 
basis {P^} such that 

P'P^ - S,, . (11) 

Like, e.g., in the context of relativity, the dual basis is 
constructed by means of some metric tensor Gij featuring 

P' = G,,P, . (12) 

A straightforward calculation yields for the inverse of the 
metric tensor 

G-j' = P^ ■ P, ■ (13) 

That is, its matrix elements are simply given by the scalar 
products of the eigenvectors of R. For the construction 
of the dual basis G^^ has still to be inverted of course. 

Within this framework we are finally able to determine 
the corrections to the eigenvalue zero. As already out- 
lined, we hereby follow precisely the scheme of standard 
time independent perturbation theory, except for replac- 
ing basis "bra" -vectors by vectors from the dual basis. 



For the eigenvalue correction (to the eigenvalue of Pqj 
i.e., zero) to first order we thus obtain 

= -iqP%Po , (14) 

where t) is a diagonal matrix in the matrix representa- 
tion used in (ITOl) with the corresponding velocities Vr, as 
diagonal elements, i.e., it is given by i?i = — iqu. If the 
system in question features full spatial mirror-symmetry 
(as our atomic wire system does) the eigenvectors Pj of 

R as well as their dual counterparts separate into two 
symmetry classes, symmetric and antisymmetric. Those 
two classes do not mix or overlap (in the sense of a regular 
scalar product). Since Pq and P° belong to the symmet- 
ric class but wPo belongs to the antisymmetric class, the 
first order correction Ei vanishes. As outlined above this 
is exactly what is required for the corresponding solution 
of the Boltzmann equation in order to describe diffusive 
dynamics at all. Thus we may indeed extract the diffu- 
sion coefficient from the second order. A straightforward 
calculation yields for the latter 

D=-^ = -P°Vr,R-^v^Po,. , (15) 

where R~^ corresponds to the inverse matrix to R ne- 
glecting the null-space. 

Thus we eventually found a diffusion coefficient de- 
scribing a diffusive class of solutions of a velocity- 
discretized, linear Boltzmann equation. Or to rephrase: 
given a specific linear Boltzmann equation we may calcu- 
late the corresponding diffusion coefficient based on (jisp . 
Thus, in our case we now have to specify a pertinent 
Boltzmann equation, i.e., the and RrjK- Of course we 
do that on the basis of the considerations in Sec. lIIIl We 
formulate a Boltzmann equation which is discretized ac- 
cording to the previously described graining. Concretely 
we identify the v^i with the slope of the dispersion re- 
lations at the centers of the corresponding grains. This 
is the standard concept of identifying particle velocities 
with group velocities. The scattering rates are taken 
from ([5]). With those allocations the Boltzmann equa- 
tion is fully specified. However, since the rates as given 
in only give rise to transitions within an energy shell, 
the dynamics completely decouple with respect to energy 
shells. Due to this decoupling the equilibrium state Pq is 
not completely determined by the Boltzmann equation, 
the latter does not fix the equilibrium particle numbers 
for full shells, it only sets the distribution onto the grains 
within one shell. Thus we have to pick a Pq in accord 
with the Boltzmann equation but additionally with rea- 
sonable overall particle numbers assigned to the shells. 
Here we do that based on the standard idea of transport 
being due to particles above the Fermi edge and due to 
holes below. Thus we weigh the equilibrium distribution 
within the shells by the probabilities of finding holes or 
particles, respectively at the corresponding energies. Do- 
ing so we find: 

Po:Po,««/(P,r)iV, (16) 
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for > and 



Po: Po^n^il- fiE,T))N, 



(17) 



for E < Sf. This construction ensures that the diffusion 
coefficient is dominated by processes occurring near the 
Fermi energy for common temperatures which appears 
physicahy reasonable. Since the diffusion coefficient as 
given in (|15|) involves the dual counterpart of the equi- 
librium distribution, an overall prefactor, i.e., the total 
number of particles does not alter its value (, hence such a 
prefactor is not specified above) . This also appears phys- 
ically reasonable. With this choice for the equilibrium 
distribution all quantities in psp are eventually speci- 
fied. 



V. NUMERICAL COMPUTATION OF THE 
DIFFUSION COEFFICIENT 

In this paragraph we concretely determine a diffu- 
sion coefficient for the model describing the atomic wire 
as introduced in Sec. by explicitly and straightfor- 
wardly evaluating formula (jl5|) . The respective compu- 
tations may be done on standard computers. The nu- 
merical effort of course depends on the finesse of the 
chosen/necessary graining, however most of the diffusion 
coefficients presented below have been found within min- 
utes of computation time. 

The calculations mainly aim at describing the typi- 
cal, characteristic dependence of the diffusion coefficient 
on two crucial system properties, i.e., the width of the 
atomic wire and the lateral coupling; we do not intend 
here to determine diffusion coefficients with extreme (nu- 
merical) precision or to describe concrete experimental 
realizations in great detail. 

Along those lines the bandwidth of the electronic dis- 
persion relation is chosen as Eb — 3eV, i.e., Ty = 0.75eV. 
Furthermore, we apply a temperature of 300K. In this 
case the equilibrium distribution Pq is mainly concen- 
trated around the Fermi energy, i.e., mainly the grains 
belonging to the shells in the vicinity of the Fermi edge 
contribute to the diffusion coefficient. Furthermore, we 
exclusively focus here on half-filling situations. For sim- 
plicity we assume that the electron-phonon coupling el- 
ements Wrs{^} from ([1]) show no strong q-dependence 
and also that there is no decisive dependence on the band 
indices r, s, i.e., we consider all momentum modes as be- 
ing equally strongly coupled, no matter which grain they 
belong to. Consequently, we set Wrs(q)^ = with 
W = O.OSeV, thus weak coupling is implemented in the 
sense of electron-phonon coupling being much smaller 
than the electronic bandwidth. Additionally, also for 
simplicity, we set the phononic dispersion to a constant 
with LOi — O.OSeV. This satisfies the above mentioned 
requirement of phonon energies being much smaller than 
electronic energies. 

As outlined above, we expect the diffusion coefficient to 
converge towards some value for finer and finer grainings. 



We clearly observe this feature in the calculations. So, 
all plotted diffusion coefficients refer to in this sense con- 
verged values, i.e., are taken from sufficiently fine grain- 
ings which here means about 100 shells per bandwidth 
Eb- This amounts to a total number of grains g « 2005. 

Fig. m shows the dependence of the diffusion coefficient 
on the lateral coupling T± for several numbers of rows (2 
to 6). 




FIG. 3: Dispersion relation for (a) 1 row, = 0.75eV, (b) 
2 rows. Til = 0.75eV, T± = 2T||, (c) 5 rows, T,, = 0.75eV, 
T, =2T,|. 



The drops of the diffusion coefficient in Fig. |4] with in- 
creasing lateral coupling may be interpreted in terms of a 
metal-semiconductor-insulator transition with respect to 
the characteristics of the underlying dispersion relation. 
Fig. 3(a) shows the dispersion relation for one row and, 
simultaneously, for any number of rows with T±_ = 0. (In 
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the latter case the band structure consists of a number 
of branches equal to the number of rows, which all coin- 
cide.) So Fi g. |3(a) refers to the utmost left data points 
in both Figs. |4] and [S] Here, the band structure is clearly 
metallic. For 2 rows the corresponding band structure is 
still metallic for Tj_ ~ T|| (see Fig. [2]), where the Fermi 
edge crosses both branches, while there is a transition to 
semiconducting behavior for Tj_ = 2T|| (see Fig. |3(b)| , 
since the Fermi energy now lies between the two bands. 
(Recall that T± determines the displacement of the differ- 
ent energy bands.) This is in accord with the decrease of 
the diffusion coefficient to very small values for T± — 2T^\ . 
Further increasing of T± would lead to the opening of a 
band gap and eventually to an insulating band structure. 

As it may be observed in Fig.|3]the diffusion coefficients 
generally decrease to small values for all depicted widths, 
which is then followed by jumps up to higher values at 
certain lateral coupling strengths, the latter varying sig- 
nificantly for different widths. The jumps occur whenever 
a band edge is driven through the Fermi edge, since the 
slope of the dispersion relation is at the band edges, 
which leads to very large rates i?,,^- This feature may be 
illustrated by the example of 5 rows with T± = 2T\\ (see 
Fig. 3(c) I, which is also covered by the data of Figs.Uand 
[51 Although the band structure is in principle metallic 
because the Fermi edge crosses one band, the diffusion 
coefficient is still very small, since the Fermi edge lies on 
the edge of two other bands, which yields high rates 
If Ti is now further increased, those two bands are soon 
driven out of the direct vicinity of the Fermi edge and 
the corresponding scattering channels (i.e., states into 
which scattering processes may occur,) get out of reach 
for not too large temperatures. This is what leads to a 
"sudden" increase of the diffusion coefficient. Especially 
note that these jumps cause a significant spreading of the 
curves above T± ss f.3r||, i.e., when the lateral coupling 
becomes stronger than the longitudinal. 

This feature is further analyzed in Fig. [SI which shows 
the resulting diffusion coefficients for different widths of 
the atomic wire, i.e., different number of parallel rows 
B, for two lateral couplings Tj_ = 2T\\ and Tj. = 
One finds that the corresponding diffusion coefficients for 
Ti = 2T|| come out rather discontinuous and become 
very small for certain widths, which may be explained 
analogously to the examples discussed above {T± = 2T\\ , 
2 and 5 rows). This effect diminishes for larger numbers 
of rows, since there are many bands which are crossed by 
the Fermi edge and only few which are tangent to it. 

For T± = T|| one finds a relatively smooth, continuous 
curve. However, in both cases we observe a dependency 
of the diffusion coefficient on the width of the atomic 
wire, which has to be classified as a finite size effect. Of 
course this effect should vanish in the limit of infinitely 
large widths B — oo, i.e., the diffusion coefficient is ex- 
pected to become constant, since the model then passes 
into a regular 2-d lattice. The plots in Fig. [5| suggest such 
a convergence at a width of approximately 50 rows and 
are therefore in accord with the above educated guess. 



Although a numerical verification was possible we omit 
such a calculation here for it would be numerically rather 
costly. 

The above argument concerning the influence of the 
properties of the electronic band structure on the inves- 
tigated dependencies of the diffusion coefficient is also 
valid for models with more realistic dispersion relations 
or electron-phonon-couplings, so that the discontinuous 
behavior will most likely remain. The same holds true for 
scattering mechanisms other than electron-phonon, such 
as weak impurity or electron-electron scattering. 




FIG. 4: Plot of the diffusion coefficient D against the lateral 
coupling T± for different widths, i.e., different numbers of 
rows, according to the legend. The plot is normalized to the 
D obtained for 1 row, denoted by Di. Other parameters: 
Eb = 3eV, Til = 0.75eV, T = 300K, haff filling. 



VI. SUMMARY AND OUTLOOK 

We investigated diffusive transport behavior of elec- 
trons in atomic wires, described on the basis of some 
pertinent quantum model. The rather simplified model is 
meant to schematically describe planar atomic wires with 
a width of a few atoms, formed, e.g., on a crystal surface 
(or embedded in some crystal bulk). The electrons within 
the wire are taken to be mutually non-interacting, diffu- 
sive behavior is solely induced by coupling to the crystal 
phonons. To this end we firstly computed a linear master 
equation describing the dynamics for some, with respect 
to momentum space, coarse grained occupation numbers, 
hereby closely following an approach from [l^. The re- 
sulting rate matrix is then interpreted as the collision 
term of a pertinent linear Boltzmann equation. Further- 
more, we analyzed solutions of linear Boltzmann equa- 
tions which can be characterized as representing diffusive 
transport dynamics. An explicit expression connecting 
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FIG. 5: Plot of the diffusion coefficient D against the number 
of rows, i.e., the width for Tx = 2^ (black dots) and T±_ = Ty 
(white triangles). The plot is normalized to the D obtained 
for 1 row, denoted by D\. Other parameters: Eb ~ 3eV, 
Ty = 0.75eV, T = 300K, half filling. 



the diffusion coefficient to the parameters of the Boltz- 
mann equation is given. Based on this formula we numer- 
ically computed diffusion coefficients for wires of different 
types. In particular, we investigated the dependency of 
the diffusion coefficient on the width of the atomic wire as 
well as the dependency on the lateral coupling strength 
between neighboring rows. We found that, if the lateral 
exceeds the longitudinal atomic coupling within the wire, 
the diffusion coefficient exhibits large jumps with respect 
to varying widths. This feature may be clearly classified 
as a finite size effect. Concretely the diffusion coefficient 
may become very small and even nearly vanish for certain 
combinations of widths and lateral couplings. There is 
also a rather strong dependence on the lateral coupling 
for fixed widths showing jumps for strong lateral cou- 
pling. The corresponding curves differ significantly sub- 
ject to the width. Thus, to summarize this results: the 
diffusion coefficient depends strongly and in a non-trivial 
way on both, the wire width and the lateral coupling. 
We additionally discussed how this discontinuous depen- 
dence may be understood in terms of the band structure 



of the non-interacting electronic system. 

By using a generalized Einstein relation [2^ the dif- 
fusion coefficients may be converted into conductivities, 
which would probably be better accessible by experi- 
ments. Therefore, the approach introduced here may 
serve as a basis for concrete transport-theoretical ana- 
lyzations to be performed in the near future. 

To obtain more realistic quantitative results, the model 
parameters used in this approach may be better adjusted 
to some concrete (experimental) realizations, e.g., a 
more realistic electron-phonon coupling, possibly includ- 
ing electron-electron interactions, impurities etc.. Addi- 
tionally, one could choose a more realistic electronic dis- 
persion relation, which may, e.g., be obtained by density 
functional theory. 

We close with the following rather speculative sugges- 
tion: One may also think of possible applications of such 
atomic wires as controlled switches or transistors on the 
nanoscale. If such an atomic wire consisting only of a 
few atomic rows is operated at a point where the Fermi 
edge and the band edge of one electronic branch coin- 
cide (as depicted, e.g., in Fig. |3(b)[ ), there is obviously 
a very strong dependence of the transport coefficient on 
small variations of the band structure, see, e.g., the case 
of B = 3, = 1.3 in Fig. H Such a variation of 

the band structure may be induced in a controlled fash- 
ion. This could be achieved by a local lateral deformation 
of the wire through some piezoelectric mechanism which 
causes a variation of the lateral coupling. More promising 
may be the concept of applying a lateral potential gra- 
dient by means of some gate electrode. Such a potential 
gradient would also have a direct effect on the relative 
positions of the electronic branches in the energy scheme 
and hence induce a strong change of the transport coeffi- 
cient. Note that such a transistor could be implemented 
without any form of doping and, just like a field effect 
transistor, would not require any gate current. 
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